The optical conductivity of graphene in the visible region of the spectrum 
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We compute the optical conductivity of graphene beyond the usual Dirac cone approximation, giving results 
that are valid in the visible region of the conductivity spectrum. The effect of next nearest neighbor hoping is also 
discussed. Using the full expression for the optical conductivity, the transmission and reflection coefficients are 
given. We find that even in the optical regime the corrections to the Dirac cone approximation are surprisingly 
small (a few percent). Our results help in the interpretation of the experimental results reported by Nair et al. 
[Science 320, 1308 (2008)]. 
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I. INTRODUCTION 

Graphene, an atomically thin material made only of car- 
bon atoms arranged in a hexagonal lattice, was isolated only 
recently. 1,2 Several reviews on the physics of graphene are al- 
ready available in the literature ! 3 ' 4 ' 5 ' 6 

At low energies, E < 1 eV, the electronic dispersion has 
the form e(fe) = ±3ifca/2, where t is the nearest neighbor 
hopping integral and a is the carbon-carbon distance. The ef- 
fective theory at these energy scales is that of a massless Dirac 
Hamiltonian in (2+1) dimensions. If the experimental probes 
excite the system within this energy range, the Dirac Hamil- 
tonian is all there is for describing the physics of graphene. 
On the other hand, for excitations out of this energy range it 
is necessary to include corrections to the Dirac Hamiltonian 
which will modify the energy spectrum and thus the density 
of states of the system. One immediate consequence is that 
the energy dispersion is no longer a function of the absolute 
value of the wave-number k. In this paper, we will calculate 
the optical conductivity of graphene including the leading cor- 
rections to the Dirac cone approximation. 

One of the first calculations of the optical conductivity of 
graphene, using the Dirac Hamiltonian were done by Gusynin 
and Sharapov. 7 This first study was subsequently revisited a 
number of times ) 8 ' 9 ' 10 and summarized in Ref. Ill 111 . How- 
ever, these authors did not include non-linear effects in the 
calculation. Also the effect of disorder was done on a phe- 
nomenological level, by broadening the delta functions into 
Lorentzians characterized by constant width T. We note that 
in the Dirac-cone approximation, the conductivity can also 
be obtained from the polarization. The calculations for finite 
chemical potential and arbitrary \q\ and ui were done by Wun- 
sch et alM and Hwang and Das SarmaJ 3 - 

The calculation of the optical conductivity of graphene, in 
the Dirac Hamiltonian limit, including the effect of disor- 
der in a self consistent way was done by Peres et al^ and 
recently also corrections due to electron-electron interaction 
were discussed. 1516 The calculation for the graphene bilayer 
with disorder was done by Koshino and Ando, 17 and by Nils- 
son et al..— The optical conductivity of a clean bilayer was 
first computed by Abergel and Falko, 19 , and recently general- 
ized to the biase d 20 ' 21 ' 22 bilayer case by Nicol and Carbotte. 23 



Within the Boltzmann approach, th e op tical conductivity of 
graphene was considered in Refs. [24,25], where the effect of 
phonons and the effect of mid-gap states were included. This 
approach, however, does not include transitions between the 
valence and the conduction band and is, therefore, restricted 
to finite doping. The voltage and the temperature dependence 
of the conductivity of graphene was considered by Vasko and 
Ryzhii, 26 using the Boltzmann approach. The same authors 
have recently computed the photoconductivity of graphene, 
including the effect of acoustic phonons M- 

The effect of temperature on the optical conductivity of 
clean graphene was considered by L. A. Falkovsky and A. 
A. VarlamoVf2i The far-infrared properties of clean graphene 
were studied in Refs. |29|and[3(l Also this study was restricted 
to the Dirac spectrum approximation. 

It is interesting to note that the conductivity of clean 
graphene, at half filling and in the limit of zero temperature, is 
given by the universal value 7re 2 /(2/i) i 31 ' 32 On the other hand, 
if the temperature is kept finite the conductivity goes to zero 
at zero frequency, but the effect of optical phonons does not 
change the value of the conductivity of clean graphene.— This 
behavior should be compared with the calculation of the DC 
conductivity of disordered graphene, which for zero chemical 
potential presents the value of 4e 2 /(7r/t)i 14 ' 31 ' 34 ' 35 

From the experimental point of view, the work of Kuz- 
menko et al. 36 studied the optical conductivity of graphite 
in the energy range [0,1] eV, and showed that its behavior 
is close to that predicted for clean graphene in that energy 
range. An explanation of this odd fact was attempted within 
the Slonczewski-McClure-Weiss model. The complex dielec- 
tric constant of graphite was studied by Pedersen for all en- 
ergy ranges. 37 The infrared spectroscopy of Landau levels in 
graphene was studied by Jiang et al^ and Deacon et al^, 
confirming the magnetic field dependence of the energy levels 
and deducing a band velocity for graphene of 1.1 x 10 6 m/s. 
Recently, the infrared conductivity of a single graphene sheet 
was obtainedi 4 ^ 4 ^ 

Recent studies of graphene multilayers grown on SiC from 
THz to visible optics showed a rather complex behavior 42 
with values of optical conductivity close to those predicted 
for graphene at infrared frequencies as well as to those mea- 
sured in graphite 36 . This experiment 42 especially indicates the 
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need for a graphene theory valid all the way to optical frequen- 
cies. The absorption spectrum of multilayer graphene in high 
magnetic fields was recently discussed in Ref. [43l including 
corrections to the Dirac cone approximation. 

In this paper we address the question of how the conduc- 
tivity of clean graphene changes when ones departs from the 
linear spectrum approach. This is an important question for 
experiments done in the visible region of the spectrum4i The 
paper is organized as follows: in Sec. [TT] we introduce our 
model and derive the current operator; in Sec. [Ill] we discuss 
the optical conductivity of graphene by taking into account its 
full density of states; in Sec. [IV] we discuss the effect on the 
optical conductivity of a next nearest neighbors hopping term; 
in Sec. [V] we analyze the scattering of light by a graphene 
plane located at the interface of two different dielectrics and 
give the transmissivity and reflectivity curves in the visible 
region of the spectrum; finally in Sec. [Vl]we give our conclu- 
sions. 



n. THE HAMILTONIAN AND THE CURRENT 
OPERATORS 

The Hamiltonian, in tight binding form, for electrons in 
graphene is written as 

H = -tJ2 E [a-UR)ba(R + S)+H.c] 

R,cr S—Si — 83 



2 ^ E [aURW(R + S)+H.c\ 

Ft, (7 S — S4 — Sq 



-|E E H(R)b,(R + S)+H.c.], (1) 

R.cr S—S4 — Sq 

where the operator a\(R) creates an electron in the carbon 
atoms of sub-lattice A, whereas b' a (R) does the same in sub- 
lattice B, t is the hopping parameter connecting first nearest 
neighbors, with a value of the order of 3 eV, and t' is the hop- 
ping parameter for second nearest neighbors, with a value of 
the order of Q.lt. The vectors Si are represented in Fig. Q~|and 
have the form 



«i = |(l,V3) , «5 2 = |(l,-V3) , 5 8 = -o(l,0) 



6 4 = a (0, Vsj , S 5 = -<5 4 , 6 6 = y f 1, 77| 



(2) 

In order to obtain the current operator we modify the hop- 
ping parameters as 

t te^ Mt) - 5 , (3) 

and the same for t' . Expanding the exponential up to second 
order in the vector potential A(t) and assuming that the elec- 
tric field is oriented along the x direction, the current operator 




FIG. 1: (color online) Representation of the vectors Si, with i — 
1 — 9. The carbon-carbon distance, a, and the A and B atoms are 
also depicted. 



is obtained from 



3a 



OH 



dA x (t) ' 

leading to j x = j£ + A x (t)j®, The operator reads 

n tie 



(4) 



J j 



E [8 x aUR)b a (R + 8)-H.c] 



R,(j S—S 1 —63 



t'ie 



-r-J2 E ISA(RW(R + S)-H.c.} 



2h 

t'ie 
~2h 



R.(T S — S4 — Sc) 



J2 E [8xbUR)b a {R+8)-H.c.].(5) 



FL,<t S—S4 — Sg 



The operator j x can be found from the linear term in A x {t) 
expansion of the Hamiltonian. 



III. THE OPTICAL CONDUCTIVITY 



A. The Kubo formula 



The Kubo formula for the conductivity is given by 



iA s (uj + i0+) ihA s (uj + i0+) ' 



(6) 



with A s = N C A C the area of the sample, and A c = 3\/3a 2 /2 
(a is the carbon-carbon distance) the area of the unit cell, from 
which it follows that 

Ra xx (u) = D6{(J) -\ — , (7) 



and 



hujA s 



<jf> ^A xx {iU + tO+) 



A s u> fkoA s 
where D is the charge stiffness which reads 



D = — TT 



<jg_> ^k xx (uj + iO+) 
A s " HA S 



(8) 



(9) 
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The function A xx (u> + i0 + ) is obtained from the Matsubara 
current-current correlation function, defined as 

rHff 

A xx {iuJ n ) = / dre l ^ T < T T f x (0) > . (10) 
Jo 

In what follows we start by neglecting the contribution of t' 
to the current operator. Its effect is analyzed later and shown 
to be negligible. The function ^sA xx (ui + i0 + ) is given by 



^-V,,U »0+) £/[#*)] 

k 



x [n F (-t\<f>(k)\ - ft) - n F (t\<f>(k)\ - /*)] 

x [tt5(uj - 2t\(j)(k)\/h) - tt(5(w + 2t\(f>(k)\/K)} , (11) 

where np(x) is the usual Fermi distribution, \i is the chemical 
potential, and the function 3tA xx (ui + i0 + ) is given by 

5RA xx (a,+i0+) = -tf^p^ fifth)] 



x [n F (-t\(/)(k)\ - (i) - n F (t\(f>(k)\ -n)] 

- (2\<t>(k)\y ' 



with 

/[</»(fc)] = 18-4|0(fc)| 2 + 18 



[jR0(fc)] 2 - [^(fc)] 
|0(fc)| a 



(13) 



and V denoting the principal part of the integral. The 
graphene energy bands are given by e(fe) = ±i|0(fe)|, with 
<j)(h) defined as 

<f>(k) = 1 + e k ^ Sl - s ^ + e fc -(<52-<5 3 ) . (14) 

B. The real part of the conductivity 

The expression for ( fTTT i can almost be written in terms of 
the energy dispersion e(fe), except for the term 



\4>(k)\ 2 



(15) 



In order to proceed analytically, and for the time being (see 
Section IlllCt . we approximate this term by its value calcu- 
lated in the Dirac cone approximation (see appendix A) 



N, 



l_ v ^(fe)] 2 -[30(fc)] 2 



\m\ 2 



g(\<fi(k)\)~0, (16) 



where g (| (f>(k) | ) is some given function depending only on the 
modulus of </>(fe). With this approximation, we have 



/[#*:)]- 18 -4|0(fc)| : 



(17) 



Introducing the density of states per spin per unit cell, p(E), 
defined as 



(18) 




u=0.2 eV, 7= 10K ' 1 

-- u=0.2eV r T=300K 






J , I , 1 i 





8 10 2 4 6 8 10 



- -r 



0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 0.5 

oj (eV) oj fsv; 



FIG. 2: (color online) The optical conductivity as function of fre- 
quency for two values of the chemical potential, /i = eV and 
fi = 0.2 eV, and two temperatures T = 10 K and T = 300 K. The 
bottom panels are a zoom in, close to zero frequency, which allow to 
depict the frequency region where differences in the chemical poten- 
tial and in temperature are most important. We have used t = 2.7eV. 



(12) the expression for the real part of the conductivity reads 



8A r huj 



p(fiw/2)[18- {hwf/t 2 



, fku + 2u Hu> — 2u 

tanh — - — — — h tanh 



4fc B T 



4fc B T 



(19) 



Equation [19] is essentially exact in the visible range of the 
spectrum; missing is only the a contribution coming from Eq. 
H~6l which contribution wil later be shown to be negligible. In 
the above equation <r is 



&0 



7T e 
2~h 



(20) 



The momentum integral in Eq. ( fT8l can be performed lead- 
ing to 



p(E) = 



2E 



where F(x) is given by 



AE/t 

y/ F (E/t)^ \F(E/t] 
V iE/t ) 



< E <t, 
t<E<3t, 



F(x) = (1 + x) 1 - 
and K(m) is defined as 



2 - I)' 



K(m) = / dx[(l ~x 2 )(l 
Jo 



-mx 2 )]- 1 ' 2 



(21) 



(22) 



(23) 



In Figure |2] we give a plot of Eq. (fT9b over a large en- 
ergy range, including the visible part of the spectrum (E £ 
[1.0,3.1] eV). 
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It is useful to derive from Eq. $1% an asymptotic expansion 
for 3ta xx (ui). For that, we expand the density of states around 
E — and obtain 



P{E) 



2E 2E 3 
+ 



10E 5 



(24) 



Using Eq. (1241) in Eq. (l9[ we obtain for the optical conduc- 
tivity the approximate result 




$t<Txx(u) = (TO 



^tanh 



72 t 2 



p + tanh - S ) • (25) 



4fc B T 



4fc B T 



In the case of jj, = this expression is the same as in Kuz- 
menko et al. 36 and in Falkovsky and Pershoguba 29 if in both 
cases the (hcu/t) 2 term is neglected. 



FIG. 3: (color online) Left: a{u))/ao as a function of the frequency, 
including both Eq. \\9\ and the correction Ka"^, Eq. {28), for sev- 
eral values of t. Right: The correction Sa" x , given by Eq. l !28t , for 
several values of t. It is clear that the contribution from this term has 
barely no effect on the results given by Eq. 1191 >. The calculations 
are for zero chemical potential and for room temperature (there is no 
visible effect on a(to)/ao in the visible range of the spectrum, when 
compared to a zero temperature calculation). 



C. Correction to $ta xx (lu) introduced by Eq. (16) 

We now want to make quantitative the effect of the term 
given by Eq. ( fT6] l, which was neglected in Eq. ( fT9l ). To that 
end we expand the function <j>(k) up to third order in momen- 
tum. The expansion is 



3a 1 / 3a 

— {k y -lk X ) + -\ — 



{kl + kl/'S + 2lk X ky) 



1 /3a 

6 It 



(ik 3 x - fc 3 /3 - 3k 2 x k y + ik 2 y k x ) . (26) 



The angular integral in Eq. ( fl6] l leads to 



7r / 3a&;\ 



where we still assume |0(fe)| = iak/2. Within this approxi- 
mation the contribution to the conductivity coming from Eq. 
( [Tol l has the form 



D. The imaginary part of the conductivity 

Neglecting the term proportional to Eq. ( fl~5b . the imaginary 
part of the conductivity is given by 

^C7 XX {U) = (7o(M--/X 3 /* logjr 7p\ 

FthJ IT 9 7T \TVjJ — ZfJ,\ 



(To f Tv~o\ . \hu) + 2/i| 
36^ \TJ ° g \huj-2fi\ ' 



(29) 



where we have included all the terms that diverge at Jko = 2/j, 
and the contribution from the cubic term in frequency in the 
density of states. The contribution of the last term of f[(f>(k)] 
in ( fT3l ) is given by 



3<» 



<jq \ ( hu> \ 



log 



18tt4!2 4 V t ) b |fej-2/i| 



(30) 



If we neglect the terms in /i 3 and oj 2 we obtain the same 
expressions as those derived by Falkovsky and Pershoguba 29 . 
We note that these terms are also obtained from the polar- 
izability in the limit q — * since the Fermi velocity is not 
/c-dependentji 2 - 



1 



= °0^4 ( — 



t J 



tanh 



fuo + 2fi 
4k B T 



tanh 



hu> — 2/i 



4fc s T 
(28) 

Due to the prefactor, this contribution has only a small effect 
and shows that the current operator basically conserves the 
circular symmetry found close to the A'-points. In Figure [3] 
we present <j(ui)/ao as function of the frequency, considering 
several values of t, in the optical range and also discussing the 
numerical value of the term given in Eq. 



E. The Drude weight and the Hall coefficient 

The Drude weight (or charge stiffness) defined by Eq. (0 
can be computed in different limits. In the case /i = we are 
interested in its temperature dependence. For zero tempera- 
ture the exact relation 



(31) 
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assures that D = when fi 
has the following form: 



D(T,fj.) = ta 



4tt 2 



\7". ^ 



0. In general, the Drude weight 



tanh 



t\<t>{k)\ + fi 
2k B T 



tanh 



8 |^(fe) 

t\m\ 



2k B T 



(32) 



In the case of finite /i, the temperature dependence of D(T, /i) 
is negligible. In the Dirac cone approximation we obtain 



D(0, fi) — AiraoH 



9 \t 



On the other hand, at zero chemical potential the temperature 
dependence of the charge stiffness is given by 



D(T,0) = 8ir\n2a k B T - 4ir((3)a 



[k B Tf 
t 2 



(34) 



where is the Riemann zeta function. 

Zotos et al. have shown a very general relation between the 
Drude weight and the Hall coefficient. 45 46 This relation is 



Rh — — "T^" 



(35) 



1 dD 

eD dn 

Equation (l35l l does not take into account the possibility of val- 
ley degeneracy and therefore it has to be multiplied by two 
when we apply it to graphene. In the case of a finite chemical 
potential we have the following relations between the Fermi 
wave vector fci? and the chemical potential: n = kp/ir and 
/j = 2takp/3. Applying Eq. ( f35T > to graphene we obtain 



R 



H 



2 n- 1 / 2 /2-3a 2 V^/8 1_ 

e \fn — a 2 n 3 / 2 /4 en 



(36) 



IV. EFFECT OF t' ON THE CONDUCTIVITY OF 
GRAPHENE 

In this section we want to discuss the effect of t' on the 
conductivity of graphene. One important question is what 
the value of t' is in graphene. Deacon et al^ proposed that 
the dispersion for graphene, obtained from a tight-binding ap- 
proach with non-orthogonal basis functions, is of the form 



E = ±- 



t\<p(k) 



'l=F«oW*)| 



(37) 



with |0(fc)| — |fca, with a the carbon-carbon distance. On 
the other hand the dispersion of graphene including t' has the 
form 



E = ±t-ka — t' 
2 



-(ka) 2 -3 
4 



(38) 



To relate t' and sq we expand Eq. (f37T > as 



E ~ ±t\<f>(k)\(l ± s \(j)(k)\) = ±t^ka + s tj(ka) 2 , (39) 



which leads t' /t — —sq with sq = 0.13. 

For computing the conductivity of graphene we need to 
know the Green's function with t' . These can be written in 
matrix form as 



G°(k,iLJ n ) = V -. 

L ^ 1.1, J„ — 



1/2 



ai|(Xfc)|/fr + 2i'[|(/>(fc)| 2 ~3]/fr 
-o0(fc)/|0(fc)| ' 



-a0(fe)7|0(fe)| 



(40) 



(33) where &°(k, iujr, ) stands for 



fpO/r. } .,\_( G A A(k,iuJ n ) G A B(k,iuJn) 
0, (K,tu n )- I G B A(k, iw n ) G BB (k,iu n ) 



(41) 



From Eq. (|40l we see that only the poles are modified, with 
the coherence factors having the same form as in the case with 
t 1 = 0. The current operator j£ — j^ t +j^ t ,, as derived from 
the tight-binding Hamiltonian is written in momentum space 
as 



Jx.t 



^]T5>(fe)-3)at >fe & CT , fe 



and 



Jx.t' 



2h 

a 

(<p*(k)- 

St'iea 



k 



(42) 



2li 



(4 



.)■ 



(43) 



The operators j£ t and j^ t , are the current operators associ- 
ated with the hopping amplitudes t and t', respectively. The 
current-current correlation function is now a sum of three 
different terms: one where we have two j£ t operators, an- 
other one where we have a and a j^ t ,, and a third 
one with two j^ t i. This last term vanishes exactly, since it 
would correspond to the current-current correlation function 
of a triangular lattice. Also the crossed term vanishes ex- 
actly which can be understood by performing a local gauge 
transformation to the fermionic operators of one sub-lattice, 
only. The first term leads to a contribution of the same form 
as in Eq. ( fT9l but with the numerators of the two tanh 
replaced by E+ = huj + 2t' [(hu) 2 / (At 2 ) - 3] + 2fi and 
E_=hw- 2t'[{huj) 2 /(U 2 ) - 3] - 2/i, respectively. 

As a consequence the effect of t' in the conductivity a 
graphene only enters in the band structure E± in the Fermi 
functions. In Figure|4]we plot the real part of the optical con- 
ductivity for two different values of \x, one with the Fermi 
energy in the conduction band and the other with the Fermi 
energy in the valence band. There is a small effect near twice 
the absolute value of the chemical potential, due to the break- 
ing of particle-hole symmetry introduced by t'. For optical 
frequencies, the effect of t' is negligible. 
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to (eV) 



FIG. 4: (color online) Real part of the conductivity for two values 
of the chemical potential at the temperature of 45 K. The parameters 
used are t — 3.1 eV and t' = — 0.13t. Only the energy range of ui £ 
[0.1, 0.3] eV is shown because only here has the chemical potential 
difference any noticeable effect. 
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FIG. 5: (color online) Geometry of p polarized light scattering be- 
tween two media with graphene separating them. The electrical per- 
mittivities of the two media are e^eo, with i = 1, 2 



V. THE ELECTROMAGNETIC SCATTERING PROBLEM 

Here we derive the reflectivity and the transmissivity of 
light between two media, characterized by electrical permit- 
tivities EiEQ, with i = 1,2, separated by a graphene flake. The 
scattering geometry is represented in Fig. [5] i.e., we assume 
the field to propagate in the direction k = (k x , 0, k z ). 

In the following, we assume the field to be given by E = 
(E x , 0, E z ) (p polarization). The case of s polarization is ad- 
dressed in AppendixlBl 

The electromagnetic boundary conditions then are 47 



(D 2 -D 1 )-n = p, 
nx(E 2 -E 1 )=0, 



(44) 
(45) 



where p is the surface charge density, in our case the graphene 
charge density. If we represent the intensity of the incident, 



reflected, and transmitted electric field as Ei, E r , and E t , re- 
spectively, the boundary conditions can be written as 

(Ei - E r ) cos 6»i = E t cos 6 2 , (46) 
-e 2 e E t sin 9 2 + eie (^ + E r ) sin 9 X = p , (47) 

where eq is the vacuum permittivity, t\ and e 2 are the relative 
permittivity of the two media and Q\ and 9 2 are the incident 
and refracted angle, respectively. Now the continuity equation 
in momentum space reads 

p(uj) = j x (u)k x /uj, (48) 

and Ohm's law is written as 

j x {u) = a(u;)E x = a(uj)E t cos6 2 . (49) 

Combining Eqs. ( |46b - j49l , we arrive at the following result, 
valid for normal incidence, for the transmissivity T 



T : 




4(e l£o ) 2 



(50) 



ei KV eie2 + £l ) e o + \/eicr(w)/c| 2 

If we now consider both media to be vacuum and that the 
graphene is at half filling (cr(w) ~ Co) we obtain 

1 



T 



(l + na/2) 2 



1 — ira , 



(51) 



where a — e 2 / (Aireoch), is the fine structure constant. The 
reflectivity is also controlled by the fine structure constant a. 
For normal incidence it reads 



R 



fela{u>)/c - eieol" 



lV e i e 2 £ + *Jeia{u)/c + eie | 
and if both media are the vacuum we obtain 



R 



-T. 



(52) 



(53) 



In Fig. |6l the transmission and reflection coefficients for 
normal incident as function of the frequency for temperature 
T = 1QK are shown where the first medium is vacuum (ei = 
1) and the second medium is either vacuum (e 2 = 1) or a 
SiC»2-substrate (e 2 — — 2, being the high-frequency 
dielectric constant of SiC^). The left hand side shows the data 
for zero doping and the right hand side for finite doping p, — 
0.2eV. In Appendix iBl we present the formulas for arbitrary 
angle of incidence. 

It is interesting to compare the result for graphene with that 
for bilayer graphene. For the bilayer, the transmissivity is 
given by 19 

T = 1 -2naf 2 (u;) 
with f 2 (o-i) given by 



/a(«) 



2(7uj + *x) {hu>/t ± ) 2 
{hhj - 2t ± )9(huj - 2t ± ) 
2{hw - tj.) ' 



(54) 
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■ \X=0, E 2 =1 
H=0, £,=2 



■ n=o, e 2 =r 

H=0, £ =2 
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— \i=0.2eV, e 2 =J 
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I 







co(el/) 



agreement between the presented theory for graphene and ex- 
periments by Dawlaty et al£^ at visible frequencies indicates 
that the interlayer interaction in epitaxial-SiC graphene is sig- 
nificant and cannot be neglected. 
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FIG. 6: (color online) The transmissivity and reflectivity for normal 
incident as function of the frequency for T = 1GK where the first 
medium is vacuum (ei = 1) and the second medium is either vacuum 
(£2 = 1) or a Si02-substrate (e2 = too — 2). Left: At zero chemical 
potential. Right: At finite chemical potential fi = 0.2eV. 



and t± the hopping amplitude between the graphene planes. 
For frequencies much larger than t±, which is the case in an 
experiment done in the visible region of the spectrum, one 
obtains 



/ 2 M~i- 



(55) 



which leads to T ~ 1 — 27ra. Again, as in graphene, the 
transmissivity is controlled by the fine structure constant. It is 
interesting to note that for Hui <C t± we also obtain the same 
result for T. 

The appearance of the fine structure constant a in the two 
cases is connected to the spinorial structure of the electronic 
wave function. In other words, the reduction of the transmis- 
sivity through a clean system is caused by a universal current 
induced by interband transitions. 



VI. CONCLUSIONS 

We have presented a detailed study of the optical proper- 
ties of graphene based on the general, non-interacting tight- 
binding model. Special emphasis was placed on going beyond 
the usual Dirac-cone approximation, i.e., we included the cu- 
bic term in the density-of-states. The conductivity was thus 
consistently calculated to order (fku/t) 2 for arbitrary chemi- 
cal potential and temperature. 

We also assessed the effect of the next nearest neighbor 
coupling t' on the optical properties. We find that the addi- 
tional terms to the current operator do not contribute to the 
conductivity and that modifications only enter through the 
modified energy dispersion. 

Using the full conductivity of clean graphene, we determine 
the transmissivity and reflectivity of light that is scattered 
from two media with different permittivity and graphene at the 
interface. Our results are important for optical experiments in 
the visible frequency range. 44 For example, the apparent dis- 



APPENDIX A: EQ. (16} UP TO FIRST ORDER IN 
MOMENTUM 

The function 4>(k) is given close to the Dirac point by 



This leads to the following result 

mk)} 2 ~ mm 2 
m - — W — 

It is now easy to see that 



(Al) 



cos(26>) (A2) 
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ddT(e)g(\<j>(k)\) = 0, 



(A3) 



where we have used the result |^(fc)| = 3afc/2, valid near the 
Dirac point. 



APPENDIX B: TRANSMISSIVITY AND REFLECTIVITY 
FOR ARBITRARY INCIDENCE 

Here we present the general formula for the transmissiv- 
ity and reflectivity of light being scattered at a plane surface 
between two media of different dielectric properties and a 
graphene sheet at the interface. 

For p polarization, the reflection and transmission ampli- 
tude are obtained from the boundary conditions of Eqs. d46i i. 
( l47j > and read 



M — 1 
M+l 



t = 



ei 2K 



e 2 M+l 



(Bl) 



with M = K + S cos 6\, where d\ denotes the incident angle 
and 



eTe c 



(B2) 



Above, k\ = ^ei^/c) 2 - k 2 x (fc* = v / e 2 (ui/c) 2 - k 2 x ) de- 
notes the perpendicular component of the incident (transmit- 
ted) wave vector relative to the interface, k x the parallel (con- 
served) component and t\ (e 2 ) is the dielectric constant of the 
first (second) medium, see Fig. [5] For s polarization, r and t 
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are independent of the angle of incident and in the Dirac cone 
approximation yield the same result as for p polarization in 
the case of normal incident (9% = 0). 

Generally, the reflection and the transmission coefficient 
are given by R = \r\ 2 and T = |i| 2 fc*/fc*, respectively. For 
a simple (non-conducting) interface, this leads to the conser- 
vation law T + R = 1. Notice that there is no such conserva- 
tion in the present case due to absorption within the graphene 
sheet. 



For a suspended graphene sheet with e% = = 1 at the 
Dirac point (<j(oj) ~ a ), the reflection and transmission co- 
efficient for p polarization read 

(gcos^) 2 1 
(l + Scos^) 2 ' (l + acos^) 2 ' K ' 

with a — ira/2 and a = e 2 / (Aireoch) the fine structure con- 
stant. 



1 K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, 
S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, Science 306, 
666 (2004). 

2 K. S. Novoselov, D. Jiang, T. Booth, V.V. Khotkevich, S. M. Mo- 
rozov, A. K. Geim, PNAS 102, 10451 (2005). 

3 A. K. Geim and K. S. Novoselov, Nature Materials 6, 183 (2007). 

4 A. H. Castro Neto, F. Guinea, and N. M. R. Peres, Phys. World 
19, 33 (2006). 

5 A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, 
and A. K. Geim, larXiv:0709.1 1631 

6 C. W. Beenakker, larXiv:0710.3848l 

7 V. P. Gusynin, S. G. Sharapov, Phys. Rev. B 73, 245411 (2006). 

8 V. P. Gusynin, S. G. Sharapov, and J. P. Carbotte, Phys.Rev.Lett. 
96, 256802 (2006). 

9 V. P. Gusynin, S. G. Sharapov, and J. P. Carbotte, Phys.Rev.Lett. 
98, 157402 (2007). 

10 V. P. Gusynin, S. G. Sharapov, and J. P. Carbotte, Phys. Rev. B 75, 
165407 (2007). 

11 V. R Gusynin, S. G. Sharapov, and J. P. Carbotte, Int. Jour, of 
Mod. Phys. B 21, 4611 (2007). 

12 B. Wunsch, T. Stauber, F. Sols, and F. Guinea, New J. Phys. 8, 
318 (2006). 

13 E. H. Hwang and S. Das Sarma, Phys. Rev. B 75, 205418 (2007). 

14 N. M. R. Peres, F. Guinea, and A. H. Castro Neto, Phys. Rev. B 
73,125411 (2006). 

15 I. F. Herbut, V. Juricic, and O. Vafek, Phys. Rev. Lett. 100, 046403 
(2008). 

16 L. Fritz, J. Schmalian, M. Mueller, and S. Sachdev, Phys. Rev. B 
78, 085416 (2008). 

17 Mikito Koshino and Tsuneya Ando, Phys. Rev. B 73, 245403 

(2006) . 

18 Johan Nilsson, A. H. Castro Neto, F. Guinea, and N. M. R. Peres, 
Phys. Rev. Lett. 97, 266801 (2006); Phys. Rev. B 78, 045405 
(2008). 

19 D. S. L. Abergel and Vladimir I. Fal'ko, Phys. Rev. B 75, 155430 

(2007) . 

20 T. Stauber, N. M. R. Peres, F. Guinea, and A. H. Castro Neto, 
Phys. Rev. B 75, 115425 (2007). 

21 Eduardo V. Castro, K. S. Novoselov, S. V. Morozov, N. M. R. 
Peres, J.M.B. Lopes dos Santos, Johan Nilsson, F. Guinea, A. K. 
Geim, and A. H. Castro Neto, Phys. Rev. Lett. 99, 216802 (2007). 

22 Jeroen B. Oostinga, Hubert B. Heersche, Xinglan Liu, Alberto F. 
Morpurgo, and Lieven M. K. Vandersypen, Nature Materials 7, 
151 (2007). 

23 E. J. Nicol and J. P. Carbotte, Phys. Rev. B 77, 155409 (2008). 

24 N. M. R. Peres, J. M. B. Lopes dos Santos, and T. Stauber, Phys. 
Rev. B 76, 073412 (2007). 



- 5 T. Stauber, N. M. R. Peres, and F. Guinea, Phys. Rev. B 76, 
205423 (2007). 

26 F. T. Vasko and V. Ryzhii, Phys. Rev. B 76, 233404 (2007). 

27 F. T. Vasko and V. Ryzhii, Phys. Rev. B 77, 195433 (2008). 

28 L. A. Falkovsky and A. A. Varlamov, Eur. Phys. J. B 56, 281 
(2007). 

29 L. A. Falkovsky and S. S. Pershoguba, Phys. Rev. B 76, 153410 

(2007) . 

30 L. A. Falkovsky, larXiv:0 806.3663 to be published in "Dubna- 
Nano2008". 

31 A. W. W. Ludwig, M. P. A. Fisher, R. Shankar, and G. Grinstein, 
Phys. Rev. B 50, 7526 (1994). 

32 N. M. R. Peres and T. Stauber, Int. J. Mod. Phys. B 22, 2529 

(2008) . 

33 T. Stauber and N. M. R. Peres, J. Phys.: Condens. Matter 20, 
055002 (2008). 

34 Nguyen Hong Shon and Tsuneya Ando, J. Phys. Soc. Jpn. 67, 
2421 (1998). 

35 Tsuneya Ando, Yisong Zheng, and Hidekatsu Suzuura J. Phys. 
Soc. Jpn. 71,1318 (2002). 

36 A. B. Kuzmenko, E. van Heumen, F. Carbone, and D. van der 
Marel, Phys. Rev. Lett. 100, 117401 (2008). 

37 Thomas G. Pedersen, Phys. Rev. B 67, 113106 (2003). 

38 Z. Jiang, E. A. Henriksen, L. C. Tung, Y.-J. Wang, M. E. 
Schwartz, M. Y. Han, P. Kim, and H. L. Stormer, Phys. Rev. Lett. 
98, 197403 (2007). 

39 R. S. Deacon, K-C. Chuang, R. J. Nicholas, K. S. Novoselov, and 
A. K. Geim, Phys. Rev. B 76, 081406(R) (2007). 

40 Z. Q. Li, E. A. Henriksen, Z. Jiang, Z. Hao, M. C. Martin, P. Kim, 
H. L. Stormer, and D. N. Basov, Nature Physics 4, 532 (2008). 

41 T. Stauber, N. M. R. Peres, and A. H. Castro Neto, Phys. Rev. B 
78, 085418 (2008). 

42 Jahan M. Dawlaty, Shriram Shivaraman, Jared Strait, Paul 
George, Mvs Chandrashekhar, Farhan Rana, Michael G. Spencer, 
Dmitry Veksler, and Yunqing Chen, arXiv:0801.3302 

43 P. Plochocka, C. Faugeras, M. Orlita, M. L. Sadowski, G. Mar- 
tinez, M. Potemski, M. O. Goerbig, J.-N. Fuchs, C. Berger, W. A. 
de Heer, Phys. Rev. Lett. 100, 087401 (2008). 

44 R. R. Nair, P. Blake, A. N. Grigorenko, K. S. Novoselov, T.J. 
Booth, T. Stauber, N. M. R. Peres, and A. K. Geim, Science 320, 
1308 (2008). 

45 X. Zotos, F. Naef, M. Long, and P. Prelovsek, Phys. Rev. Lett. 85, 
377 (2000). 

46 P. Prelovsek and X. Zotos, Phys. Rev. B 64, 235 1 14 (2001). 

47 John David Jackson, Classical Electrodynamics, (Wiley, 2001, 
3 rd edition), p. 16. 



